Structure of a superconducting vortex pinned by a screw dislocation 
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Spatial dependence of the magnetic field and the superconducting current in a flux line pinned by a 
screw dislocation are computed. Interaction of a superconducting vortex with the chiral-symmetry 
breaking elastic strain of a screw dislocation results in a helical current along the axis of the dislo- 
cation. It is argued that screw dislocations make impossible a force-free arrangement of flux lines in 
the presence of a transport current. 

PACS numbers: 74.20. De, 74.60.-w 
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Pinning of vortices by screw dislocations can be of in- 
terest because of th.&,|experimental evidence of this ef- 
fect in Y-Ba-Cu-OBBua Studies of pinning by disloca- 
tions are almost as old as the Abrikosov's idea of su- 
perconducting Mortices (see, e.g., the review by Camp- 
bell and EvettsB and references therein). More recently, 
these studies have been revived in connection with high- 
temperature superconductors They are normally based 
upon the expansion of scalar parameters of a supercon- 
ductor, such as, e.g., Tc, in terms of the deformation 
tensor and its trace. 

A superconducting vortex coupled to a screw dislo- 
cation is an interesting theoretical problem, because of 
the broken chiral symmetry of the deformation field pro- 
duced by such a dislocation. Same as the magnetic field 
B, a screw dislocation running, e.g., along the Z-axis, 
is an axial vector, say p, that is invariant under reflec- 
tion. Consequently, p = p-z = ±1 is a pseudoscalar 
that changes sign under reflection. This makes possible 
a linear relation between the polar vector of the current 
density j and the axial vector of the magnetic field B: 
jocpB, the relation that would be otherwise prohibited 
by the invariance with respect to reflections. 

A closely related problem of a vortex coupled to-a spiral 
defect has been studied by Ivlev and ThompsonQ'B for an 
extreme case of a layered superconductor with a Joseph- 
son coupling between the layers. In such spiral 
defect, running perpendicular to the layers, geometrically 
connects them by a continuous helical path around the 
defect. Ivlev and Thompson elegantly solved this prob- 
lem in spiral coordinates and, in accordance with the 
above symmetry arguments, demonstrated the existence 
of a "fountain-like" current along the axis of the defect. 

The purpose of this paper is to solve the problem in the 
opposite extreme case of a flux line coupled to a screw dis- 
location in an isotropic three-dimensional superconduc- 
tor. In such case the existence of longitudinal currents 
parallel to the dislocation core is somewhat less obvious. 
The model we suggest is complementary to the model of 
Ivlev and Thompson. It is based upon the description 
of dislocations within continuous elastic theory. We find 



that superconducting currents do flow along screw dislo- 
cations in three-dimensional superconductors, although 
the spatial dependence of these currents is different from 
the one found in Ref.7. Our results should be directly 
relevant to MgB2 and conventional isotropic supercon- 
ductors. The Y-Ba-Cu-0, where the coupling between 
vortices and screw dislocations has been experimentally 
observed, falls somewhere between the two models. 

The free energy of the system is — Td + J^n + J^gl, 
where Td is the energy of the dislocation, Tn is the 
energy of the normal electron liquid in the absence of 
the magnetic field, and Tql is the Ginzburg-Landau free 
energy. 
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Here A is the vector potential (B = VxA), i/j = 
I'i/'l exp(i</)) is the complex order parameter of the su- 
perconducting phase, a and b are constants, Di is the 
gauge-invariant derivative. 
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and TOy is the tensor of effective masses. For an isotropic 
superconductor, in the absence of crystal defects, rriik — 

The presence of a dislocation results in a non-zero elas- 
tic strain. 
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At distances exceeding a few lattice spacings from the dis- 
location core- the components of the dimensionless tensor 
are smalla and the parameters of the superconductor, 
such as a, 6, and m^, can be expanded into the power 
series of Uij. We shall see that for a screw dislocation 
Tr(uij) = 0, that is, screw deformations change the sym- 
metry of the crystal but not the local density. Thus, to 
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the lowest order in Uij, the interaction of the screw dis- 
location with the Ginzburg-Landau order parameter can 
be introduced by the substitution 
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where g is a dimensionless parameter. From the physics 
of the electron states in crystals, it is clear that the ef- 
fect of the lattice deformations on the tensor of effective 
masses is small as long as Uij are small. If some of Uij 
are large, the effect should be also large. Consequently, 
g must be of order unity. 

Since crystal defects are insensitive to superconductiv- 
ity, Td must have a very weak dependence on ip. We will 
neglect that dependence. Then the variation of Tgl with 
respect to and A gives the following two Ginzburg- 
Landau equations: 
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and the Maxwell equation V x B 
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where gjj = 5ij + gu^ . 



Eq. (g) describes the vortex core where the concentra- 
tion of Cooper pairs, li/jp, changes from zero at the center 
of the vortex to a constant value, |a|/6, at distances ex- 
ceeding the coherence length ^ = h/2{m\a\Y^'^ . At such 
distances the second Ginzburg-Landau equation becomes 



Ai5,^i(VxB), 



A, 



2tt 



(7) 



where (j) is the phase of V', $o = hc/2e is the flux quan- 
tum, and \l — (mc^/STre^j-fAp)^^^ is the London penetra- 
tion depth. We shall consider the case of Al >> ^ which 
is relevant to high-temperature superconductors. At a 
large distance from the vortex core, where B exponen- 
tially goes to zero, Eq. reduces to A = (<i>o/27r)V0. 
After the integration over a closed distant contour enclos- 
ing the vortex, it produces the conventional condition of 
the quantization of the magnetic flux. This condition 
remains unchanged by the deformations. Applying curl 
to both sides of Eq. , one obtains a modified London 
equation 



X^^B, + eyfcVj-5j:;i(VxB), = 0. 
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We can now turn to the problem of a flux line centered 
at a screw dislocation parallel to the Z-axis of the crystal. 
It should be naturally studied in circular cylindrical co- 
ordinates: z, r, Lf. At distances, r, exceeding a few lattice 
spacings from the core of the screw dislocation, the only 
non-zero component of the displacement field u iaj 
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where p = ±1 is the chirality of the screw dislocation 
and b is the Burgers vector that coincides with the lattice 
spacing in the Z-direction. Consequently, the only non- 
zero components of the elastic strain are 
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Although the linear elastic theory fails near the axis of 
the dislocation, for the purpose of our study the above 
formulas are exact as long as the coherence length ^ and 
the London penetration depth Al are greater than 6. As 
is shown below, the significant effect of the dislocation 
on the superconducting vortex comes from the fact that 
the clastic strain given by Eq. ( |lo|) is long range. 

It is convenient to introduce a dimensionless coordi- 
nate p = r/Ai and a dimensionless small parameter 
(3 = gb/inXL- Substituting Eq. @ into Eq. (|), and 
retaining terms of the lowest order in /3, one obtains the 
following two equations for non-zero components of the 
magnetic field Bz{p) and B^{p): 
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For a flux line parallel to Z, the condition /? << 1 results 
in B^ « Bz for all r > ^. Consequently, the effect of 
the dislocation on Eq.(ll) can be neglected and the sys- 
tem of equations for Bz{p) and B^{p) can be solved by 
iteration. Putting /3 = in Eq.(ll) one obtains a con- 
ventional solution for the Z-component of the magnetic 
field inside the flux line. 
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where Kq is the modified Bessel function and the coeffi- 
cient in front of Kq is such that the total magnetic flux 
through the XY-plane equals <I>o. Then Eq.(12) can be 
reduced to the following form 
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The solution of this equation, that goes to zero at r— >oo, 
is 
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where C is a constant of integration (to be computed 
later). 

The current density in the flux line is given by j = 
(c/47r)VxB. The (/^-component of j is the conventional 
vortex current 
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FIG. 1. Radial dependence of the longitudinal current den- 
sity in a superconducting vortex coupled to a screw disloca- 
tion. 
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The unusual feature of the problem, as in Ref.7, is the 
presence of the Z-component of the current, 
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Substituting here of Eq. mm , we obtain 
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The longitudinal current jz changes sign at r of order A^. 
At r >> Al it becomes exponentially small. For a screw 
dislocation of the opposite chirality, the current is in the 
opposite direction. 

A screw dislocation either ends at the surface of the 
crystal or forms a loop inside the crystal. In the first 
case, the boundary condition, UigijDj^ — (with n 
being the unit vector normal to the boundary), prohibits 
currents through the surface of the superconductor. This 
results in a surface current that flows outward from the 
vortex core at ^ < r < Al. In the case of a dislocation 
loop, the flux line will follow the loop. 

The coupling of the flux line to a screw dislocation 
is the (6/C)^ fraction of the vortex-core energy avail- 
able for pinning. In high-temperature superconductors, 
where the coherence length can be of order of the Burg- 
ers parameter, dislocations can provide a r|athe*; strong 
pinningcl, in accordance with observationsu'0'LlQ. One 
should keep in mind, however, that the presence of the 
B^p component of the field in the flux line pinned by a 
screw dislocation makes impossible the "force-free" situ- 
ation in which the field is parallel to the transport cur- 
rent. In that sense columnar pins with no chirality have 
the advantage over screw dislocations. 

I thank Lev Bulaevskii and Tolya Kuklov for useful 
discussions. This work has been supported by the U.S. 
Department of Energy through Grant No. DE-FG02- 
93ER45487. 
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Note that the dependence of on p, given by Eq. ([TS 
differs from the one obtained by Ivlev and Thompson! 

The magnetic field at distances significantly exceeding 
Al is not possible in a superconductor. Consequently, 
the total superconducting current J^, that flows parallel 
to the dislocation outside the vortex core, must be zero. 
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where we have introduced k = A^/C- With the help of 
Eq. (|l^) and Eq. (|l8|), this condition allows one to obtain 
constant C: 



C = ln(K) + kKq{1/k)K^^{\/k). 
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A few observations are in order. The n dependence 
of iz{p) is rather weak: C changes from C = 0.1916 
at K = l/'\/2 (which is the boundary of type-II super- 
conductivity) to the asymptotic form: C — 21n(K), at 
K^oo. For K = 100 and p = 1, jz{f) is shown in Fig. 
|l|. The maximal longitudinal current, jz ~ {b/£,)jo (with 
jo = c<I'o/12-\/37r2f A| being the Ginzburg-Landau criti- 
cal current) occurs at r ~ ^ near the core of the vortex. 
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